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We present a tight binding theory of the Dirac point resonances due to adsorbed atoms and 
molecules on an infinite 2D graphene sheet based on the standard tight binding model of the graphene 
7r-band electronic structure and the extended Hiickel model of the adsorbate and nearby graphene 
carbon atoms. The relaxed atomic geometries of the adsorbates and graphene are calculated using 
density functional theory. Our model includes the effects of the local rehybridization of the graphene 
from the sp 2 to sp 3 electronic structure that occurs when adsorbed atoms or molecules bond cova- 
lently to the graphene. Unlike in previous tight-binding models of Dirac point resonances, adsorbed 
species with multiple extended molecular orbitals and bonding to more than one graphene carbon 
atom are treated. More accurate and more general analytic expressions for the Green's function 
matrix elements that enter the T-matrix theory of Dirac point resonances than have been available 
previously are obtained. We study H, F, OH and O adsorbates on graphene and for each we find 
a strong scattering resonance (two resonances for O) near the Dirac point of graphene, by far the 
strongest and closest to the Dirac point being the resonance for H. We extract a minimal set of tight 
binding parameters that can be used to model resonant electron scattering and electron transport 
in graphene and graphene nanostructures with adsorbed H, F, OH and O accurately and efficiently. 
We also compare our results for the properties of Dirac point resonances due to adsorbates on 
graphene with those obtained by others using density functional theory-based electronic structure 
calculations, and discuss their relative merits. We then present calculations of electronic quantum 
transport in graphene nanoribbons with these adsorbed species. Our transport calculations capture 
the physics of the scattering resonances that are induced in the graphene ribbons near the Dirac 
point by the presence of the adsorbates. We find the Dirac point resonances to play a dominant 
role in quantum transport in ribbons with adsorbates: Even at low adsorbate concentrations the 
conductance of the ribbon is strongly suppressed and a transport gap develops for electron Fermi 
energies near the resonance. The transport gap is centered very near the Dirac point energy of 
for H, below it for F and OH and above it for O. We find conduction in ribbons with adsorbed H 
atoms to be very similar to that in ribbons with equal concentrations of carbon atom vacancies. We 
predict ribbons with adsorbed H, F, OH and O, under appropriate conditions, to exhibit quantized 
conductance steps of equal height similar to those that have been observed by Lin et al. [Phys. Rev. 
B 78, 161409(R) (2008)] at moderately low temperatures, even for ribbons with conductances a few 
orders of magnitude smaller than 2e 2 /h. 
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I. INTRODUCTION 

In recent years graphene nanoribbons have been 
the subject of increasing experimental-^- and 
theoretical^— interest. Ideal ribbons that are 
uniform in width and free of defects, adsorbates and 
other disorder should transmit electrons ballistically 
(i.e. without scattering) and consequently, as is the 
case for other quasi-one-dimensional (ID) ballistic 
nanostructures^ their low temperature conductances 
are expected to be quantized in integer multiples of 
2e 2 /h^^^L-^r^l The ribbons that have been 

realized experimentally to date have been far from 
ballistic; consistent with theoretical work on strongly 
disordered ribbon3^ r 27 ' 33 i 38 ' 40 i 40 " — their measured 
conductances when a few transverse subbands are 
populated with electrons have been much smaller than 
the conductance quantum 2e 2 //i. Thus the recent exper- 
imental observation of quantized conductance steps in 
graphene nanoribbons by Lin et al* 3 - at moderately low 



temperatures was surprising and especially so in view of 
the fact that the observed conductance step heights were 
two orders of magnitude smaller than 2e 2 /h. In a pre- 
vious paper j 4 * we explained this puzzling phenomenon 
as arising from enhanced electron backscattering near 
subband edge energies due to the presence of defects. 
This explanation is consistent with the conclusion 
drawn by Lin et alj£ that the conductance steps that 
they observed were evidence of subband formation 
in their nanoribbon samples. In the models that we 
studied 4 ^ carbon atom vacancies in the interior of the 
ribbon played a crucial role: They were responsible 
for the formation of equally spaced conductance steps 
in a range of temperatures T high enough to suppress 
universal conductance fluctuations but for which ksT is 
smaller than the subband spacing, in agreement with the 
experiment. 3 Here fc# is Boltzmann's constant. However, 
whether such vacancies were actually present in the 
experimental samples at the required concentrations was 
not determined in the experimental work of Lin et ali^ 
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and it is widely believed based on STM 51 and TEM 52 
measurements that graphene samples can be free from 
carbon atom vacancies over large areas. On the other 
hand it is reasonable to expect atomic and molecular 
species to be adsorbed on graphene ribbons prepared 
using presently available fabrication techniques ! 41 ' 53 " — 
Thus it is of interest to explore the possible role that 
such adsorbates may play in conductance quantization 
of the kind reported by Lin et al4 

A remarkable property of pristine graphene is that 
the electron dispersion near the Fermi energy is linear 
and forms Dirac-like cones in fc-space centered on two 
points K and K' in the Brillouin zone^ The energy 
at which the density of states of graphene vanishes is 
known as the Dirac point. It has been suggested that im- 
purities that strongly perturb the graphene should give 
rise to resonant states in the vicinity of the Dirac point 
and that these resonant states result in strong scattering 
of electrons in the graphene^r— Adsorbed atoms and 
molecules that are covalently bonded to graphene per- 
turb the graphene strongly and thus may be expected to 
give rise to such "Dirac point resonances" . However, how 
adsorbate-induced Dirac point resonances affect electron 
transport in graphene nanoribbons is a topic that is yet 
to be explored theoretically or experimentally. 

In this paper we generalize the previously proposed 
analytic theories of the impurity-induced Dirac point 
resonances^— to the case of adsorbates on graphene 
whose electronic structure is described within the ex- 
tended Hiickel model of quantum chemistry^ The ex- 
tended Huckel models is a semi-empirical tight bind- 
ing scheme that provides a simple but reasonably re- 
alistic description of the electronic structures of many 
molecules. It has been used successfully in explaining 
and predicting experimental transport properties of a 
variety of molecular systems 50 including conduction in 
molecular wires bridging metal contacts 7 ^— and molec- 
ular arrays on silicon^ 9 - - — and the electroluminescence,— 
current-voltage characteristics 8 ^ and STM images^ 3 , of 
molecules on complex substrates. Thus it offers a nat- 
ural way to extend the standard tight binding model of 
pristine graphene^ 9 - to the case of graphene with adsor- 
bates. 

Here we develop a tight binding model of graphene 
with adsorbates based on extended Huckel theory and 
use it to carry out quantum transport calculations for 
graphene nanoribbons with adsorbates. An important 
advantage of our approach is that it includes multiple 
atomic or molecular orbitals of the adsorbed species as 
well as the relevant non-7r graphene orbitals explicitly in 
the tight binding model used to calculate the effect of 
Dirac point resonances on transport whereas other tight- 
binding transport calculation o 64 i 70 ' 71 (which have been 
for 2D graphene with adsorbates or impurities) have been 
restricted to much simpler smg/e-orbital models of Dirac 
point resonances. 

As specific examples of adsorbates we consider H, F 
and O atoms and OH groups, species that may have 



been present in the experimental samples of Lin et al^ 
that were made by oxygen plasma reactive ion etching 
using a hydrogen silsesquioxane etch mask that was later 
removed in a hydrofluoric acid solution. 

We estimate the relaxed geometries of these adsorbates 
and carbon atoms to which they bond using ab initio den- 
sity functional calculations^ 4 - It is energetically favorable 
for the graphene carbon atoms to which adsorbed species 
bond to move out of the graphene plane towards the ad- 
sorbate by fractions of an Angstrom and partial rehy- 
bridization of the carbon atom from sp 2 to sp 3 bonding 
then occurs i 41 ' 53 i 68 ' 86 Therefore in our treatment of the 
Dirac point resonances we include the atomic valence or- 
bitals of the carbon atoms that are involved in the sp 3 
bonding in addition to the atomic valence orbitals of the 
adsorbed species and the 2p z orbitals of the graphene car- 
bon atoms that are included 5 ^ in the standard tight bind- 
ing model of graphene. Our theory also applies to species 
that bond simultaneously to more than one graphene car- 
bon atom as is the case for an adsorbed O atom. 

We develop more accurate and more general analytic 
expressions than have been available to date^Sr— for 
the matrix elements of the Green's function of pristine 
graphene that enter the theory of the Dirac point res- 
onances and check their accuracy by numerical calcula- 
tions. We then calculate the matrix elements of the T- 
matrix that describes the scattering of graphene electrons 
by the adsorbate as a function of energy and thus deter- 
mine the energies at which the Dirac point resonances 
of the various adsorbates occur as well as the resonance 
energy profiles. 

For each of the H, F and OH adsorbed species we find 
a strong resonance located near the Dirac point. For the 
adsorbed O atom we find the T-matrix to exhibit a more 
complex energy profile with a pair of overlapping reso- 
nances of different widths near the Dirac point. For each 
of these adsorbed species we also develop a minimal set of 
tight binding parameters that yield an accurate descrip- 
tion of its Dirac point resonance(s). These parameter 
sets are used in our transport calculations on graphene 
nanoribbons that we report here and are expected also to 
be useful in other theoretical work such as studies of the 
tunneling spectra of adsorbates on graphene that may 
be observed in scanning tunneling spectroscopy experi- 
ments. 

We show that Dirac point resonances due to adsor- 
bates have a strong signature in the transport charac- 
teristics of graphene nanoribbons that depends strongly 
on the adsorbed species even at low adsorbate concentra- 
tions. We investigate the possible influence of adsorbates 
on the quantized conductances that have been observed 
experimentally 3 in graphene nanoribbons at moderately 
low temperatures, and that we have studied theoretically 
previously^ - in models of graphene nanoribbons with va- 
cancies. The results obtained here regarding adsorbate- 
induced Dirac point resonances in graphene also provide 
a physical interpretation of key features of the calculated 
nanoribbon conductances that we shall present. The 
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transport calculations that we present yield the following 
salient results:: 

(i) The adsorption of each of the species that we study 
on a graphene ribbon results in very strong electron scat- 
tering especially at energies in the vicinity of the Dirac 
point scattering resonance associated with the respective 
adsorbate. This in turn leads to strong suppression of 
the ribbon conductance and a transport gap opening up 
in a range of electron Fermi energies near the energy of 
the Dirac point resonance. 

(ii) The transport gaps occur for electron Fermi ener- 
gies around the Dirac point for H, below the Dirac point 
for F and OH and above the Dirac point for O adsor- 
bates. The case of O differs qualitatively from those of 
F and OH because O binds to two carbon atoms belong- 
ing to different graphene sublattices while the F and OH 
bind to a single carbon atom. 

(iii) The conductance characteristics as a function of 
the electron Fermi energy for ribbons with adsorbed H 
atoms are very similar both qualitatively and quantita- 
tively to those of ribbons with the same concentration of 
carbon atom vacancies. 

(iv) We predict that ribbons with each of these ad- 
sorbed species should under appropriate conditions ex- 
hibit equally spaced conductance steps at moderately low 
temperatures even for adsorbate concentrations for which 
the conductance is much smaller than 2e 2 /h, consistent 
with experiment £ 

The remainder of this article is organized as follows. 
In Sections [TTI- IVIII we present a theory of the infinite 2D 
graphene sheet with adsorbed atoms or molecules. Then 
in Sections I Villi and |IX] we apply the results obtained 
in the preceding sections to the problem of transport in 
graphene nanoribbons. We formulate our tight binding 
model of the electronic structure of adsorbed atoms and 
molecules on graphene in Section [II] In Section Mil we 
show how the tight-binding Hamiltonians can be trans- 
formed into effective graphene Hamiltonians that include 
the effect of the adsorbate by generalizing previous the- 
oretical work to the case of adsorbate species with mul- 
tiple molecular orbitals and bonding to more than one 
graphene carbon atom. In Section IIVI we briefly discuss 
the T-matrix theory used to study Dirac point resonances 
in graphene analytically. In Scction|V]we derive an exact 
relation between the graphene Green's function matrix 
elements that enter the T-matrix theory. We use this 
relation to obtain accurate analytic expressions for these 
matrix elements and check their accuracy by exact nu- 
merical calculations. In Section lVll we present our results 
for the Dirac point resonances of H, F and O atoms and 
OH groups on graphene as well as a minimal tight bind- 
ing model that accurately reproduces those results and 
is used in the transport calculations that we present in 
Section HXl In Section lVll we also examine quantitatively 
the effect of the adsorbate-induced sp 3 rehybridization of 
the graphene on the Dirac point resonances and find it 
to be important, and especially so for adsorbed H atoms. 
In Section IVIII we compare our results for the properties 



of Dirac point resonances due to adsorbates on graphene 
with those obtained by others using density functional 
theory-based electronic structure calculations. In Sec- 
tion |VHl] we discuss the model and methodology used in 
our calculations of transport in graphene nanoribbons. In 
Section llXl we present: In Section llX Al we discuss ribbons 
with H adatoms including the effects of Dirac point reso- 
nances and graphene rehybridization on transport and on 
the electronic density of states. Ribbons with F, OH and 
O adsorbates are then considered in Sections IIXBI IIX CI 
and I IX PI We comment on conductance asymmetries of 
ribbons relative to the Dirac point that result from the 
presence of adsorbates in Section IIX El In Section IIX Fl 
we focus on the topic of conductance quantization in rib- 
bons with adsorbates. We discuss renormalization of the 
energies of the ribbon subbands and Dirac point due to 
the presence of adsorbates in Section IIX Gl and the de- 
pendence of the transport gaps on the adsorbate concen- 
tration in Section llX HI In Section [X] we summarize our 
main findings and comment briefly on potential relevant 
experiments. 



II. THE TIGHT BINDING HAMILTONIAN 

Our starting point is the simplest tight-binding model 
of pristine graphene 5 - embodied in the Hamiltonian 

Here —t is the Hamiltonian matrix element between 
nearest-neighbor 2p z carbon orbitals of the graphene lat- 
tice and a\ is the creation operator for an electron in 
2p z carbon orbital i. This Hamiltonian with t = 2.7 eV 
is known to describe the 7r band dispersion of graphene 
well at energies around the Dirac point ) 59 i 88 i.e., in the 
energy range of interest in the present work. We extend 
this tight binding model to include the adsorbate and its 
coupling to the graphene carbon atoms in the following 
way. 

We performed ab initio geometry relaxations based on 
density functional theory for the adsorbed species on 
the honeycomb graphene lattice using the Gaussian 09 
software packaged The relaxed geometries calculated 
in this way are expected to be accurate since density 
functional theory has been well optimized for carrying 
out accurate ground state total energy calculations on 
which these relaxations are based^ The structures stud- 
ied were graphene disks of several tens of carbon atoms 
passivated at the edges with hydrogen, the adsorbed 
species being bonded to the graphene near the center 
of the disk. The atoms of the adsorbed species and the 
carbon atoms with which they bond were allowed to re- 
lax freely, the other carbon atoms being held fixed in 
the standard hexagonal graphene geometry with the C- 
C distance of 1.42 A. The relaxed structures obtained in 
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FIG. 1: (color online) Relaxed geometries of adsorbates on 
graphene. C, H, F, and O atoms are black, blue, green and 
red respectively.— (a)Adsorbed hydorgen atom. H atom and 
C atom to which H binds are 1.47 and 0.35 A above graphene 
plane, (b) Adsorbed fluorine. F atom and C atom to which F 
binds are 1.83 and 0.36 A above graphene plane, (c) Adsorbed 
hydroxyl group. H atom, O atom and C atom to which O 
binds are 2.78, 1.83 and 0.41 A above graphene plane, (d) 
Adsorbed oxygen. O atom and C atoms to which O binds are 
1.51 and 0.27 A above graphene plane. The C atoms to which 
the O binds are separated by 1.47 A. The x-axis is parallel to 
the line joining the C atoms to which the O binds. 

this way are shown in Fig. [T}£^ The tight binding model 
Hamiltonian Eq. ([1]) was extended to include the atomic 
valence orbitals of the adsorbate and their coupling to the 
valence orbitals of the graphene carbon atoms by calcu- 
lating the relevant matrix elements within the extended 
Hiickel models 

Extended Hiickel theory is formulated in terms of 
small basis sets of Slater- type atomic orbitals {|</>i)}, 
their overlaps $y = (4>i\4>j), and a Hamiltonian matrix 
Hij — (<Pi\H\^j)- The diagonal Hamiltonian elements 
Ha = £i are chosen to be the experimentally determined 
atomic orbital ionization energies £j. In the present work 
the nondiagonal elements are approximated as in Ref. Iz2 
by Hij = (1.75 + Afj - O.r.-.A.UV, , (£i + £j) /2, where 
Aij = (£i — £j)/(£i+£j), a form chosen^ for consistency 
with experimental molecular electronic structure data. 
In the standard tight binding Hamiltonian ([1} of pristine 
graphene, the energy scale is chosen so that the carbon 
2p z orbital energy is zero whereas in extended Hiickel 
theory^ the carbon 2p z orbital energy is the ionization 
energy £q , = —11.4 eV. Accordingly, for consistency, in 
our extended Hiickel Hamiltonian matrix we make the 
replacement Ha — > Ha — £c Pz ■ Because the extended 
Hiickel basis states on different atoms are not in general 
mutually orthogonal the non-diagonal extended Hiickel 
Hamiltonian matrix elements are then also adjusted ac- 
cording to 

Hij -> Hij - S i: j£c Pz (2) 

as is discussed in Refill. 

Let Hf, and be the extended Hiickel Hamiltonian 
and overlap matrices defined in this way but restricted to 
the Hilbert subspace R spanned by the valence orbitals of 
the adsorbate and valence orbitals of the graphene other 
than the 2p z graphene orbitals that in the present model 
are already described by Hq. In the numerical results 
presented in this paper in addition to the adsorbate va- 
lence orbitals we include in R the 2s, 2p x and 2p y orbitals 
(see FigO} of the carbon atom(s) to which the adsorbed 



atom or molecule bonds and of its three nearest carbon 
atom neighbors. We calculate Hf, and Sj}; for the relaxed 
geometries shown in Fig. [T] We then solve the extended 
Hiickel Schrodinger equation 

H K 4> a = e a S K i> a (3) 

(for a single adsorbed H, F or O atom or OH group) nu- 
merically for its eigenstates ip a and energy eigenvalues 
e a . The eigenstates ip a obtained in this way are mutu- 
ally orthogonal. They should be regarded as extended 
molecular orbitals (EMOs) of the adsorbate. They are 
linear combinations of the atomic valence orbitals of the 
adsorbate and some of the atomic valence orbitals of the 
graphene (as is detailed above Eq. [3]) but do not include 
any graphene 2p z atomic orbitals. The EMOs will play 
a central role in the theory that follows. 

The EMOs ip a together with the 2p z orbitals of the 
graphene carbon atoms form the basis set for our tight 
binding Hamiltonian H of the graphene-adsorbate sys- 
tem that we write in the form 

H = H + ^2 t a did a + locj (dlaj + h.c.) , (4) 

Here aj is the destruction operator for an electron in the 
2p z orbital <fij of carbon atom j. d) a is the creation oper- 
ator for an electron in EMO ip a that is an eigenstate of 
Eq. ([3]) and e a is the corresponding energy eigenvalue. 
laj — (ipalHlcfij) the matrix element of the extended 
Hiickel Hamiltonian between the 2p z orbital <f>j of car- 
bon atom j and EMO ip a . For simplicity, in this paper 
we include in the sum over j in the last summation on the 
RHS of Eq. (0]) only the 2p z orbital of the carbon atom 
that is closest to the adsorbed moiety (or in the case of 
the adsorbed O the closest two carbon atoms), and we 
also neglect any changes in t in Eq. |T]) that occur due 
to the change in the graphene geometry induced by the 
adsorbate. The latter latter effect is however taken into 
account in the numerical nanoribbon transport calcula- 
tions that we report in Section [IXI 

We note that with the above definitions, the couplings 
between all of the valence orbitals (2s, 2p x , 2p y and 2p z ) 
of the carbon atom to which the adsorbate bonds and all 
of the 2s, 2p x , 2p y and 2p z valence orbitals of that carbon 
atom's nearest carbon atom neighbors are included in the 
present model that is summarized by Eq. In this 

way we include in our calculation all of the carbon atom 
valence orbitals that participate in the adsorbate induced 
sp 3 bonding. 

Because the basis set used in extended Hiickel theory 
is non-orthogonal the overlap a a j = (il> a \(f)j) between the 
2p z orbital <f>j of carbon atom j and EMO ip a may be 
non-zero. This overlap is neglected in Eq.(|3]). It has 
been show n 90 ' 91 that transport problems formulated in a 
nonorthogonal basis can be solved by transforming to an 
alternate Hilbert space in which the basis is orthogonal 
but the effective Hamiltonian matrix elements become 
energy-dependent. This transformation is the foundation 
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of the standard methods used today to treat basis set 
non-orthogonality throughout the molecular electronics 
transport literature. For the present system, the trans- 
formation is accomplished by replacing j a j in Eq. ^ by 
laj — £&aj- Here e is the electron energy at which the 
Landaucr electron transmission probability through the 
system is calculated. This correction is included in the 
numerical results that we present in this paper for the 
adsorbate induced Dirac point resonances, although in 
the interests of clarity it will not appear explicitly in the 
formulae that we present in the remainder of this article. 

III. EFFECTIVE HAMILTONIANS 

In the simplest possible model of an adsorbate repre- 
sented by just one atomic orbital a that couples only 
to the 2p z orbital of only one carbon atom j of the 
graphene the tight binding Hamiltonian of the graphene 
and adsorbate is Hi = Hq + ead^da + ^ a j (<i Q aj + h.c.) 
where the notation is as in Eq. 0] The eigenstate 1^) 
of H\ with energy eigenvalue e can then be written as 
|*) = I*,) + |*o) where |* g ) and |* a ) are the projec- 
tions of onto the space spanned by the 2p z orbitals of 
graphene and onto the orbital of adsorbed atom, respec- 
tively. With these definitions, it has been shown^i that 
l^g) is an exact eigenstate of an effective Hamiltonian 
if ff = iio + VjOjCij with the same energy eigenvalue e 

as Here Vj = J^j/i 6 ~ e «)- Thus for the purpose 

of calculating the transport coefficients of graphene with 
such an adsorbed atom within Landauer theory, it is suf- 
ficient to replace the Hamiltonian Hi with H e s, i.e, the 
Hamiltonian of graphene without the adsorbed atom but 
with an energy dependent potential 7„,-/ (e — e a ) on car- 
bon atom j of the graphene sheet. 

We note that a general theory of systems with one or 
more discrete states coupled to a continuum of states was 
presented by Fano in 1961. 9? Graphene with an adsorbed 
atom or molecule is such a system. Fano's analysis of 
such systems^ starts in the same way as the analysis in 
Ref. [6J that we have outlined above by projecting the 
eigenstates of the system onto the continuum and discrete 
state manifolds. However, unlike in Ref. [64], Fano did not 
reformulate the problem in terms of an effective Hamil- 
tonian H e ff that acts on the continuum subspace only^ 
Subsequently, effective Hamiltonians have been employed 
to study bound states coupled to continua but those ef- 
fective Hamiltonians have been non-Hermitian operators 
obtained by eliminating the continuum subspace from the 
theory^ unlike the Hermitian effective Hamiltonians ob- 
tained in Ref. [64] (and in the theory presented below) by 
eliminating the discrete state subspace. 

In the present work we need to include more than 
one extended molecular orbital ipa per adsorbed moi- 
ety in the tight binding Hamiltonian given by Eq. Q 
for the adsorbed H, F or O atom or OH group. This is 
the case even for H (which has only one valence orbital 
in extended Hiickel theory) because we include several 



graphene atomic orbitals in the subspace R in which we 
calculate the EMOs tp a for the H adsorbate as is discussed 
111 Section HU 

We find that the argument presented in Ref. that 
leads to the effective Hamiltonian H c g discussed above 
can be generalized in a direct way to adsorbates for which 
more than one effective molecular orbital and/or bonding 
of the adsorbed moiety to more than one carbon atom (as 
in the case of adsorbed O) must be considered. 

For adsorbates with more than one extended molecu- 
lar orbital that bond strongly to a single carbon atom, 
we find that the effective Hamiltonian still has the form 
H c g = Ho + Vja^aj but the effective potential Vj becomes 

Vj = J2 a l7cy | 2 /( e ~~ e a) where the sum is over the ex- 
tended molecular orbitals a of the adsorbed moiety that 
bonds to carbon atom j of the graphene. 

For a single adsorbed O atom that bonds to two neigh- 
boring graphene carbon atoms 1 and 2 the effective 
Hamiltonian is H e g = Ha + Vna\ai+V22a\a2 + Vi2aia2 + 
V 2 ia\ai where V nm = J2 a 7cm 7am/ ( e _ £ a) and the sum- 
mation is over the extended molecular orbitals a of the 
O adsorbate. For the adsorbed O atom we consider 22 
EMOs in this paper. They are linear combinations of 
the O 2s, 2p x , 2p y and 2p z valence orbitals and the 2s, 
2p x and 2p y valence orbitals of each of the six carbon 
atoms shown in Fig. [Hd) ■ The effect of the adsorbed O 
atom on the Hamiltonian eigenstates projected onto the 
graphene ir band subspace is equivalent to the combined 
effect of energy dependent potentials applied to the 2p z 
valence orbitals of the two carbon atoms to which the 
O atom bonds and an energy dependent change in the 
Hamiltonian matrix elements between those 2p z carbon 
orbitals. 

The preceding results for the effective Hamiltonians 
H e g apply equally to adsorbates on graphene nanorib- 
bons or on 2D graphene. We now use them to develop a 
better understanding of how H, F, OH and O adsorbates 
resonantly scatter electrons by extending the general T- 
matrix approach considered previously in Refs l60ll62| - 
[6^l6^l67ll69ll7ll . We consider here for simplicity the case 
of an isolated single H, F, OH or O atom or molecule on 
infinite 2D graphene. In Section IIXI we will relate our 
findings to the results of our numerical transport calcu- 
lations for graphene nanoribbons. 

IV. T-MATRIX FORMALISM 

The T-matrices that we consider are defined in the 
standard way by 

G = G° + G°TG° (5) 

where G = (e + i-q — H^)^ 1 is the full Green's func- 
tion based on the effective Hamiltonians H c s discussed 
above for a single adsorbed atom or molecule, G° = 
(e + ir] — Ho) is the unperturbed Green's function for 
7r band electrons in clean graphene and T characterizes 
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the scattering strength due to the adsorbate. T can be 
written in the standard form 



T = V + VG°V + VG°VG°V + 



(6) 



where V = a^aj J2 a \laj\ 2 /( e — O for a H, F, OH atom 
or molecule with EMOs a bound to carbon atom j. For 
an O atom with EMOs a bound to two neighboring C 
atoms 1 and 2, V = Viia\ai + V22a\a2+Vi2a\a2+V2\a\ai 

where V nm = J2 a l an l* am l ( e ~ £ a)- 

Taking matrix elements of Eq. [6]between the graphene 
2p z orbitals of the carbon atom(s) to which the adsorbed 
atom or molecule binds and summing the resulting series 
yields 



f = ( 1 - VG° 



V 



(7) 



where for the O atom adsorbate T, V, G° and 1 are the 
2x2 matrices (m\T\n), (m\V\n), (m\G°\n), and S mn with 
\m) and \n) being the 2p z orbitals of the carbon atoms 
to and n to which the O atom bonds. Here to = 1,2 
and n = 1,2. For H, F and OH that bond to one C 
atom (labelled 1) f, V, G° and 1 are the scalars (1|T|1), 
(1|V|1), (1|G°|1), and 1 respectively. 



V. ANALYTIC EXPRESSIONS FOR THE 
MATRIX ELEMENTS OF G° 

The unperturbed Green's function matrix elements 
(m\G°\ri) = G° m between the relevant 2p z carbon or- 
bitals of 2D graphene that enter Eq. ([7]) are given by 



r° (A = 

^ rn n \ C J / _j 



k,p 



(m\<I> kp )(<S> kp \n) 
e + irj- ($ kp \H \<f> kp ) 



(8) 



where \§ kp ) are the eigenstates of the unperturbed 
Hamiltonian Hq of 2D graphene given by Eq.([T|), with 
wave vector k and band index p = ±1. 

We evaluate the diagonal matrix element (1|G°|1) = 
G^e) as follows: For m = n = 1 and small e the sum- 
mand in Eq. [3] is strongly peaked in fc-space around 
the Dirac points. We therefore approximate the sum 
over the Brillouin zone in Eq. [5] by the sum of inte- 
grals over two circles in fc-space centered on the two 
Dirac points K and K', choosing the area of each cir- 
cle to be equal to half of that of the Brillouin zone. We 
then evaluate the integrals by linearizing ( < I > fcp|i?o| < I > fcp) 
in fc about each Dirac point in the standard way which 
yields ($> kp \Ho\Q kp ) s» ±3ir|fc|/2 where r is the nearest 
neighbor spacing between graphene carbon atoms. This 
yields^ 



G? x (e) 



-h 



Wit 2 



(9) 



For comparison we carried out an exact numerical eval- 
uation of Eq. [5] without linearizing (<& kp \Ho\$ kp ) or ap- 
proximating the unperturbed model Hamiltonian Hq = 



— J2(n, m ) t { a \i a m + h.c.) of 2D graphene or its eigenval- 
ues or eigenfunctions in any other way. We found the 
analytic approximation (Eq. ^ to agree with our exact 
numerical results to within a few percent for small values 
of e. We also found that the accuracy of Eq. d9]) can be 
improved and extended to larger |e| with the help of em- 
pirical correction factors a(e) = 1.07(1 + 0.66e 2 /£ 2 ) and 
/3(e) = l + 0.31e 2 /i 2 + 0.33e 4 /i 4 . The resulting expression 



G?i(e) 



ea(e) 

v^Trf 2 "'" \V3nt 2 ~e 2 



In 



V3£ 2 



(10) 



is accurate in the range \e\/t < 0.8. We note that sev- 
eral less accurate analytic approximations for G^ (e) have 
also been proposed in the literaturei 60 i 63 ' 66 i 69 ' 71 The first 
these^ underestimated G^e) by more than a factor of 
2, due in part to the use of an inadequate model of the 
graphene electronic structure with only a single Dirac 
cone instead of the two such cones that are centered at 
the Dirac points K and K' of graphene in reality. The 
most accurate of the m 60 i 63 i 66 i 69 i 71 is that given by Eq. 
37a in Ref. with e replaced by |e| in the imaginary 
term. G^ n (e) can also be expressed exactly in terms 
of elliptic integrals that, however, must be evaluated 
numerically^ 

The analytic approximation scheme that leads to Eq. 
injis inappropriate for evaluating G°x(e) directly because 
in that case the summand in Eq. (|8l) is not maximal at 
the Dirac points K and K 1 where the linear approxima- 
tion to (& kp \Ho\$ kp ) is accurate. However an analytic 
expression for G^e) that is accurate in the same energy 
range as Eq. (IT0|) can still be obtained as follows. 

We start from the identity 



(l\(e + ir)-H )G Q (e)\l) = 1 
from which it follows that 



(11) 



(e + ir ? )<l|G (e)|l) + ^<l|(-ff )|n><n|Go(e)|l) = 1 

n=2 

(12) 

In Eq. [TT1and[T2l |1) is the 2p z orbital of a carbon atom 
and |n) for n = 2,3,4 are the 2p z orbitals of its three 
nearest neighbors. Analysis of Eq. [8] shows that be- 
cause of the three-fold rotational symmetry of infinite 2D 
graphene (2|G (e)|l) = (3|G (e)|l) = (4|G (e)|l). Then, 
since (1|(— H^)\n) = t and taking the limit rj — » we 
obtain from Eq. [12] the exact result that 

G° 1 (e) = (2|G (e)|l) = |-^G? 1 (e) (13) 

It is then straight forward to show also that G^e) — 
G? 2 (e)- Finally, inserting Eq. [TUlinto Eq. [13] we obtain 



G 2 i(e) = ^7 —in 



(14) 



3t 3V3tt* 3 V\/37rf 2 -e 2 / 3V3t 

We have compared the analytic expression (fl"4"|) with the 
results of our exact numerical evaluation of Gt^e) and 
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found Eq. (fT4|) to be accurate under the same conditions 
as Eq. (fTU|) . as expected. Previous theoretical wor k 62 i 66 
has only yielded analytic expressions for G^e) for the 
special case e = 0. For that case Eq. reduces to 

G?i(0) = ji which agrees with the result stated in Ref. 
l66l but differs in sign from that stated in Ref. l62j It 
should be noted that the sign of G^ is not arbitrary; 
the calculated Dirac point resonance for O adsorbed on 
graphene is modified significantly if this sign is reversed. 



VI. DIRAC POINT RESONANCES OF H, F, OH 
AND O ADSORBATES ON GRAPHENE 



10 5 


I I 


I I 


10 4 


H 




10 3 






10 2 


OH /V 
f >x / J 


\Ao 


10 1 






10° 


I I 


i i v 



-0.75 -0.5 -0.25 0.25 0.5 0.75 



£ 

FIG. 2: (color online) Calculated square modulus of the T- 
matrix vs. electron energy e for an H, F or O atom or OH 
group adsorbed on graphene in the geometries shown in Fig. 
ff] T and e are in units of t = 2.7eV. The Dirac point of 
graphene is at e = 0. For H, F and OH T = (1|T|1). For 
O the square of the Frobenius norm of the matrix (m\T\n) 
is plotted. The EMOs included in this calculation are linear 
combinations of the atomic valence orbitals of the adsorbed 
species and the 2s, 2p x and 2p y valence orbitals of each of 
the carbon atoms shown in Fig. [T]for the respective adsorbed 
species. Thus the local rehybridization of the graphene horn 
the sp 2 to sp 3 bonding is included in the model. The overlaps 
o a j between the EMOs and the 2p z orbitals of the carbon 
atoms to which the adsorbed moieties bond are included in 
the calculations. 

The strength of scattering associated with a defect is in 
general proportional to the square modulus of appropri- 
ate matrix elements of the T-matrix. Thus the energies 
e at which resonant scattering by H, F, OH and O ad- 
sorbates should occur are those at which |(m|T|n)| 2 have 
maxima. We find these energies to be close to those at 



which |1 - VG°|- 2 for H, F and OH or |dei(l - VG )|~ 2 
for O have maxima. 

The square moduli of the matrix elements of the T- 
matrix defined in Section IIVI and calculated using the 
tight binding parameters ~f a j and e a obtained from ex- 
tended Huckel theory as is described in Section |H] are 
shown vs. the electron energy e in Fig. [5] The molecular 
orbitals ip a included in these calculations are linear com- 
binations of the atomic valence orbitals of the adsorbed 
species and the 2s, 2p x and 2p y valence orbitals of all 
of the carbon atoms shown in Fig. [T] for the respective 
adsorbed species. 13, 16, 17 and 22 EMOs are included 
in the calculations for H, F, OH and O respectively. The 
overlaps a a j between the EMOs ip a and the 2p z orbitals 
<f)j of carbon atoms j to which the adsorbed moieties 
bond are included in the calculations as is discussed at 
the end of Section [II] 

For each adsorbed species the T-matrix displays a 
prominent resonant peak (a double peak for O in Fig. 
[2]) in the vicinity of e = 0, the Dirac point of graphene. 
The electron energy edr at which the resonance is 
centered depends on the adsorbed species. £dr = 
-0.136t, -0.089*, -0.0026* for F, OH and H, respectively. 
For O there is a narrow peak near 0.112* that overlaps a 
broader peak centered near 0.090*. The strengths of the 
F, OH and broader O resonances (as measured by the 
area under the resonance curve when plotted on a lin- 
ear scale) are all comparable. The narrow O resonance 
is an order of magnitude weaker than these while the H 
resonance is two orders of magnitude stronger. 

For comparison the results of a similar calculation 
but for a simpler model in which the carbon atom 
2s, 2p x and 2p y valence orbitals are omitted from the 
EMOs are shown in Fig. [3] In this case edr = 
-0.222*, -0.194*, -0.138* and 0.046* for F, OH, H and 
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FIG. 3: (color online) Calculated square modulus of the T- 
matrix vs. electron energy e for a simpler model of H, F or 
O atom or OH group adsorbed on graphene than that used 
in the calculations presented in Fig. [2] Here only the va- 
lence orbitals of the adsorbed species themselves (no carbon 
orbitals) are included in the EMO's. The parameters e a and 
y a j used are from Table ITR the molecular orbitals ip a are used 
in the case of OH. The overlaps <r a j are also included in the 
calculation. Notation as in Fig. O 
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TABLE I: Minimal set of effective tight-binding parameters 
e a and7 aj in units of t = 2.7 eV for adsorbed H, F and and 
OH on graphene. The EMO energies t a are measured from 
the Dirac point energy of graphene. ± means that 7^ has 
opposite signs for the two carbon atoms to which the O atom 
bonds. 



adsorbate 








H 


-0.0383 


2.219 




F 


-10.862 
-2.460 
-0.914 


4.363 
1.645 
1.180 




OH 


-8.536 
-1.820 
-0.709 


3.203 
1.779 
1.540 




O 


-5.356 
-1.448 
-0.373 
0.130 
1.463 


3.240 
±1.000 

1.095 
±0.176 

1.650 





O, respectively^ Comparing Fig. [5] with Fig. [3] it is ev- 
ident that the coupling of the adsorbate to the graphene 
carbon 2s, 2p x and 2p y valence orbitals (that are in- 
volved in the partial rehybridization of the carbon atom 
to which the adsorbate bonds to the sp 3 electronic struc- 
ture) can affect adsorbate- induced Dirac point resonances 
very strongly: It is directly responsible for the H reso- 
nance in Fig. [2] being two orders of magnitude stronger 
than the resonances for F, O and OH. It is also respon- 
sible for the double peak structure of the O resonance in 
Fig. [5] that is absent in Fig. [3] Notice also the antireso- 
nance in the O T-matrix near e = 0.55i in Fig. [5] that is 
absent in Fig. [3l 

The Dirac point resonance energy increases from F to 
OH to H to O in both models and the sign of the res- 
onance energy for each species (negative for F, OH and 
H and positive for O) is the same in both models. How- 
ever, the coupling of the H adsorbate to the graphene 
carbon 2s, 2p x and 2p y valence orbitals results in the H 
Dirac point resonance being extremely close to the Dirac 
point of graphene (cdr = — 0.0026<) in Fig. [2] As will 
be seen in Section [IXI this results in the conductances of 
graphene nanoribbons with hydrogen adsorbates being 
almost symmetric about the graphene Dirac point (as 
they are for ribbons with carbon atom vacancies^), in 
marked contrast to the asymmetric conductances for the 
other adsorbed species. 

Finally we find that an accurate description of the 
Dirac resonance profiles in Fig. [5] (including their en- 
ergies, widths and heights) can be obtained by includ- 
ing in the tight binding Hamiltonian Eq. (|4]) relatively 
small sets of EMOs since some of the EMOs couple only 
weakly to the graphene ir system. The EMO energies e a 
and coupling parameters 7^ for a minimal tight binding 



model Hamiltonian that describes the Dirac point reso- 
nance profiles shown in Fig. [5] are presented in Table HI 
The EMO energy and coupling parameter sets presented 
in Table fl] will be used in the more sophisticated version 
of our transport calculations on graphene nanoribbons 
that we report in Section IIXI The values of e a given 
in Table U are the energies of EMOs calculated from ex- 
tended Hiickel theory as discussed in Section [TTJ However 
the values of some of the j a j that are given have been 
adjusted so that the small set of parameters e a and j a j 
given in Table U yields a good fit to the resonance pro- 
files in the energy range —0.75t < e < 0.75t shown in 
Fig. [2] without the need to include the overlaps o~ a j in 
the calculation. 



VII. COMPARISON WITH DENSITY 
FUNCTIONAL THEORY-BASED MODELS OF 
DIRAC POINT RESONANCES 

In the previous attempts to construct tight binding 
models of Dirac point resonances for use in transport 
calculations in graphene, the tight-binding parameters 
were obtained by fitting the results of ab initio density 
functional theory-based electronic structure calculations 
to very simple tight binding models. In those models the 
adsorbed atom or molecule M was described by just a 
single effective orbital energy parameter cm and a single 
coupling parameter 7m- For example, en = 0.66i and 
7h = 0.22t were obtained for adsorbed hydrogen and 
£oh = — 2.9t and 70H = 2.3t were obtained for the ad- 
sorbed hydroxyl group in Ref. Ifril On the other hand, 
in Ref. l7Cj much smaller values of the effective orbital 
energy parameter cm < O.lt were found for a number of 
covalently bonded adsorbed species, including hydrogen, 
and 7m > 2t was found for the same species. 

That such very different results have been obtained 
from density functional theory-based calculations even 
for hydrogen, the simplest of all adsorbed species, raises 
the question whether density functional calculations, al- 
though putatively a "first principles" method, are a 
sound basis for theoretical studies of the Dirac point res- 
onances of graphene with adsorbates. This demonstrated 
lack of consistency may be related to the fact that density 
functional theory, although well suited to calculations of 
the total ground state energies of many condensed matter 
systems, is known to have important fundamental defi- 
ciencies as a methodology for electronic structure calcu- 
lations; for a recent discussion of the relevant physics and 
a review of the literature the reader is referred to Ref. 

As a consequence of these deficiences, for example, den- 
sity functional calculations underestimate the band gap 
of silicon and other semiconductors by as much as a factor 
of two. They also yield offsets between the energy levels 
of molecules adsorbed on silicon and the silicon valence 
band edge that more sophisticated theories indicate to be 
in error by as much as lAeVMZ> The latter error is similar 
in size to the above-mentioned discrepancy between the 
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values of en predicted by the density functional theory 
calculations reported in Refs. [64| and [lO for H adsorbed 
on graphene. 

Because different density functional theory-based elec- 
tronic structure calculations have yielded such different 
results even for the Dirac point resonance due to hydro- 
gen on graphene, we chose instead to base our electronic 
structure calculations on the well known semi-empirical 
extended Hiickel model of quantum chemistry. The ex- 
tended Hiickel model has been parameterized 7 -^ based on 
a large body of experimental electronic structure data for 
atoms and molecules and has been used successfully in 
electronic transport calculations for a variety of molecu- 
lar systems, as has been outlined in Section fl] The rel- 
ative merits of electronic structure calculations based on 
density functional theory and those based on extended 
Hiickel theory have been discussed in detail in Section 
4.7 of Rcf. HO- 

The minimal tight-binding model parameters for the 
Dirac point resonance due to hydrogen adsorbed on 
graphene that we derived from extended Hiickel theory 
using the methodology described in Sections HIllVIl are 
given in Table HI They are in excellent agreement with 
the corresponding density functional theory-based results 
£m < O.li and 7m > 2t stated in Ref. 70: but are not con- 
sistent with the density functional theory-based results 
e H = 0.66i and 7 H = 0.22t reported in Ref. [H 

We find one-orbital tight binding models such as those 
used in Refs. H3 and [70 not to yield a satisfactory de- 
scription of Dirac point resonances for the other adsorbed 
species (F, OH and O) that we have considered within 
extended Hiickel theory. Therefore comparing our tight- 
binding models for the Dirac point resonances of those 
species directly with the corresponding tight binding pa- 
rameters that have been derived from density functional 
theory (as we have done above for the case of adsorbed 
hydrogen) is not possible. However, we have calculated 
the local densities of states (LDOS) associated with the 
Dirac point resonances within our extended Hiickel based 
model for these adsorbates and compare them below with 
the corresponding LDOS features calculated using den- 
sity functional theory in Ref. l68l . 

Our calculated local densities of states for graphene 
with an adsorbed H, F or O atom or OH group are shown 
in Fig. HI The extended Hiickel theory-based model 
used in these calculations is the same as that used to 
calculate the square moduli of the T-matrices that are 
shown in Fig. O That is, the EMOs included the cal- 
culations are linear combinations of the atomic valence 
orbitals of the adsorbed species and the 2s, 2p x and 2p y 
valence orbitals of each of the carbon atoms shown in 
Fig. Q] for the respective adsorbed species. Thus the 
model includes the effects of the local rehybridization of 
the graphene due to the presence of the adsorbate. The 
quantities plotted in Fig. [4] are the partial LDOS defined 
by D n (e) — — Im((n|G(e)|n))/7r where G is defined by 
Eq. [5] and |n) represents a 2p z orbital of a carbon atom. 
The dashed black curves are for graphene carbon atoms 



to which the adsorbed moieties bond and the solid red 
curves are for nearest (carbon) neighbors of those car- 
bon atoms. Note the different energy scale used for the 
case of H and also that the dashed black curve for H has 
been scaled up for clarity by a factor of 500. In each plot 
the Dirac point is at e = 0. The energies at which the 
peaks of the densities of states occur in Fig. @] are rea- 
sonably close to the energies of the Dirac point resonance 
peaks for the respective adsorbates in Fig. [2j the largest 
discrepancy being ~ 20% for fluorine. 

The partial LDOS in Fig. @]for H, F and O can be com- 
pared with the corresponding density functional theory- 
based results in Figs. 1(c), 1(d) and 5(a) of Ref. |68|, 
respectively. However, some differences between the sys- 
tems considered should be noted: The results in Fig. 0] 
are for a single adsorbed atom or molecule on an infinite 
graphene sheet whereas those in Ref. [H are for periodic 
structures with 4x4 graphene supercells each containing 
an adsorbed moiety. Thus in Ref. l68l . in contrast to the 
present work, the Dirac point resonance is expected to 
be broadened due to the presence of multiple adsorbed 
atoms on the graphene. Also, the DOS features for H, 
F and O in Ref. [68| are located relative to the Fermi 
energy. The latter may be close to the Dirac point en- 
ergy but its location relative to the Dirac point energy 
is not determined precisely since in the model systems 
studied in Ref. [68] there are no very large regions of pris- 
tine graphene with no adsorbate where the Dirac point 
is well defined. 

With these caveats, the LDOS for the Dirac point res- 
onance for H in Fig. @] is consistent with that in Ref. 
[H, although for the latter the resonance is broader and 
the partial LDOS on the carbon atom to which the H 
bonds is not as weak relative to the partial LDOS on 
its nearest carbon atom neighbors. The LDOS for the 
fluorine Dirac point resonance in Ref. |68| is also similar 
to that in Fig. 4 although the LDOS peaks in the latter 
are somewhat lower in energy relative to the Dirac point 
than those in the former are relative to the Fermi energy. 
However, as we have already noted, the precise location 
of the Fermi level relative to the Dirac point is uncertain 
in the density functional theory-based calculations. 

The LDOS for O in Fig. [4] differs markedly from that 
in Ref. |68f There is no obvious peak in the LDOS in the 
immediate vicinity of the Fermi energy in the latter case 
in contrast to the peaks associated with the Dirac point 
resonance ~ O.lt above the Dirac point in the former. 
The reasons for this difference are not clear at present. 
We note, however, that, as has been discussed above, dif- 
ferent density functional theory based calculation o 64 i 70 
have yielded very different results even for the Dirac 
point resonance for adsorbed hydrogen. Furthermore, 
recent theoretical works^ has shown even gross features 
of the electronic structures of narrow graphene nanorib- 
bons with high concentrations of adsorbed O at the rib- 
bon edges (including the presence or absence of a large 
band gap at the Fermi level) calculated using density 
functional theory to be sensitive to the precise choice of 
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FIG. 4: (Color online) Calculated local densities of states (LDOS) vs. electron energy e for an H, F or O atom or OH group 
adsorbed on graphene in the geometries shown in Fig. [1] The LDOS is in units of 1/t and e is in units of t = 2.7eV. For F, 
O atom OH the black dashed curves show the LDOS associated with the 2p z orbital of a carbon atom to which the adsorbed 
moiety bonds. For H the black dashed curve shows the LDOS associated with the 2p z orbital of a carbon atom to which the 
adsorbed atom bonds multiplied by 500. The red solid curves show the LDOS associated with the 2p z orbital of a carbon atom 
that is a nearest neighbor of a carbon atom to which the adsorbed moiety bonds. The Dirac point of graphene is at e = 0. The 
EMOs included in this calculation are linear combinations of the atomic valence orbitals of the adsorbed species and the 2s, 
2p x and 2p y valence orbitals of each of the carbon atoms shown in Fig. Q] for the respective adsorbed species. Thus the local 
rehybridization of the graphene from the sp 2 to sp 3 bonding is included in the model. 



the exchange-correlation energy functional used in the 
calculations. While the calculated LDOS for O in Ref. 
|68| resembles qualitatively that found for a "double impu- 
rity" in a simple tight binding model 68 it should be noted 
that in that model the "double impurity" represents two 
atoms that are adsorbed on adjacent carbon atoms and 
do not interact with each other directly, a situation that 
is very different than a single oxygen atom that bonds to 
two adjacent carbon atoms in its lowest energy configura- 
tion in reality and in our extended Huckel theory-based 
model. 

Given the differing theoretical predictions that we have 
discussed above, it is evident that experiments probing 
adsorbate-induced Dirac point resonances in graphene 
would be of considerable interest. 



VIII. MODEL HAMILTONIAN FOR 
TRANSPORT CALCULATIONS IN GRAPHENE 
NANORIBBONS 

In this Section we describe how the tight-binding 
Hamiltonians developed in the preceding Sections for ad- 
sorbates on infinite 2D graphene are adapted for calcu- 
lations of quantum transport in graphene nanoribbons 
with adsorbates that we consider in the remainder of this 
paper. 

We describe the graphene ribbons by the tight-binding 
Hamiltonian 

H = H 1X + Y^ £ a did a + laj (4 a i + h.c) (15) 
where 

H v = — Uj \ a\cij + h.c^j . (16) 

(hj) 



tij is the Hamiltonian matrix element between nearest- 
neighbor 2p z carbon orbitals of the ir band of the 
graphene nanoribbon. a\ is the creation operator for an 
electron in 2p z carbon orbital i. S a creates an electron 
in an extended molecular orbital (EMO) ip a that is as- 
sociated with an adsorbed atom or molecule and has en- 
ergy e a . As defined in Section HH an EMO is a linear 
combination of the valence orbitals of the adsorbed atom 
or molecule and (in the more sophisticated versions of 
the model) the 2s, 2p x and 2p y valence orbitals of the 
graphene carbon atom(s) to which the adsorbate bonds 
and neighboring graphene carbon atoms. The graphene 
2p z orbitals are not included in the EMOs since they are 
included in the tight binding Hamiltonian H through H K . 
In the present work we will include for simplicity only 
coupling matrix elements jaj between the EMOs asso- 
ciated with the adorbates and the 2p z valence oribitals 
of the graphene carbon atoms to which that adsorbed 
atoms or molecules bond. 

In the graphene ribbon 7r-band Hamiltonian we in- 
clude nearest neighbor Hamiltonian matrix elements Uj- 
For most of these we set Uj = t = 2.7eV, the usual value 
for tight binding theories of pristine graphene^iS, 8 How- 
ever, interaction with the H, F, OH and O adsorbates 
shifts the carbon atoms to which these moieties bond 
out of the graphene plane by fractions of an Angstrom, 
perturbing the values of Uj between the carbon atoms 
to which the adsorbates bond and their neighbors. Al- 
though we find the effect of this change in Uj on electron 
transport in ribbons to be modest (typically less that a 
15% difference in the conductance) we include it in the 
calculations presented in this article, estimating the mod- 
ified values of Uj by applying extended Huckel theory to 
the relaxed geometries of the graphene in the presence 
of the adsorbates £2- We do not, however, consider edge 
reconstruction effects and spin and electron interaction 



11 



phenomena; these are outside of the scope of the present 
study. 

An important effect associated with the shifting of car- 
bon atoms out of the graphene plane that occurs upon 
adsorption of H, F, OH or O is the local partial rehy- 
bridization of the graphene from sp 2 to sp 3 bonding. We 
shall elucidate the role that the rehybridization plays in 
transport in graphene ribbons by comparing the results 
of transport calculations for two models: 

(i) A simpler model in which the EMOs associated with 
adsorbed atoms or molecules are approximated by linear 
combinations of only the valence orbitals of the adsorbed 
species. 

(ii) A model in which the EMOs are linear combina- 
tions of the the valence orbitals of the adsorbed species 
and of the appropriate valence orbitals (2s, 2p x and 2p y ) 
of the graphene carbon atoms involved in the rehybridiza- 
tion, i.e., the carbon atoms to which the adsorbate bonds 
and their nearest neighbors. 

For case (ii) the number of EMOs per adsorbed moiety 
is too large for exact quantum transport calculations to 
be carried out with the computational resources available 
to us for ribbons of experimentally relevant sizes and ad- 
sorbate concentrations of interest. However it was found 
in Section lVll that reduced sets of between 1 and 5 EMOs 
per adsorbed moiety with suitably adjusted values of "/ a j 
are sufficient to provide an accurate description of the 
Dirac point resonances in graphene induced by H, F, OH 
and O adsorbates. These minimal sets (with parameters 
listed in Table HJ) will be used in the transport calcula- 
tions for case (ii) that we present in this paper. 

In our transport calculations the adsorbed atoms and 
molecules are introduced by randomly placing them on 
the graphene surface. They are characterized by the 
probability p to find an adsorbed moiety per carbon 
atom. In order to convert p to the usual concentration 
one should scale it by the number of carbon atoms di- 
vided by the sample size, 3.8 x 10 19 m~ 2 . 

In the linear response regime the zero temperature con- 
ductance of the graphene ribbon is given by the Landauer 
formulaSSLSa-iSa 

G=^J>. (17) 

ij 

where Tji is the transmission coefficient from the subband 
i in the left lead to the subband j the right lead, at the 
Fermi energy. For non-zero temperatures 

where T(E) = YUj Tji(E) and f(E) is the Fermi distri- 
bution function. Tji is calculated by the recursive Green's 
function method, see Ref. [H for details. 



TABLE II: Tight-binding parameters for the adsorbed H, F, 
O and OH in units of t — 2.7 eV for the simplest model in 
which the EMOs do not include the carbon 2s, 2p x and 2p y 
atomic orbitals. The atomic and molecular orbital energies 
e a are measured from the Dirac point of graphene. 7a,c 2pz 
is the Hamiltonian matrix element between the orbital ip a of 
the adsorbate and the 2p z orbital of the nearest graphene C 
atom. The ± means that the 'y a ,c 2pz values for the two C 
atoms to which the O atom bonds have opposite signs. For 
OH the parameters are given for both the atomic O and H 
orbitals and the molecular orbitals ij) a of OH. 7c 2pz ,c 2pz is 
the Hamiltonian matrix element between the 2p z orbital of 
the graphene C atom to which the adsorbed moiety bonds 
and the 2p z orbital of a neighboring C atom.— For O the first 
7c 2pz ,c 2pz is for the pair of C atoms to which the O bonds 
while the second is for a C atom to which the O bonds and 
another nearest C neighbor of that C atom. 



adsorbate 


1p<x 






lc 2vr ,c 2vr 


H 


Is 


-0.81 


1.89 


0.79 


F 


2s 
2p z 


-10.59 
-2.48 


4.70 
1.45 


0.79 




2s° 


-7.74 


4.10 




OH 


2p? 


-1.26 


1.24 


0.73 
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IX. RESULTS 

In our transport calculations the ribbon has a width 
W = 30nm as in the experiments of Lin et ah— The ad- 
sorbed H, F, OH or O are assumed to be present only in 
a hnite region of length L — 500nm of the ribbon which 
is attached at its two ends to semi-infinite leads repre- 
sented by ideal ribbons of the same width. The edge 
configuration is taken as armchair in the following. For 
comparison we also present some results for ribbons with 
interior carbon atom vacancy defects4^ 



A. Ribbons with H Adatoms 

1. Atomic Geometry and Electronic Structure 

Hydrogen is the simplest adsorbate that bonds cova- 
lently to a carbon atom of the graphene lattice. Ac- 
cording to the extended Hiickel model, the H Is orbital 
energy locates not far from the Dirac point of graphene, 
ejj ls = —0.814, and the Hamiltonian matrix element be- 
tween the hydrogen Is orbital and the carbon 2p z orbital 
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FIG. 5: (color online) Calculated conductances as a function of the Fermi energy for graphene ribbons with different adsorbed 
species at a concentration p — 10~ 4 . The relaxed geometries of the adsorbed species and nearby graphene atoms are shown 
in the insets. The dotted grey lines show the conductance of the ideal ribbon without any defects. The solid black lines in 
(a), (c), (e) and (g) show the calculated conductances of ribbons with H, F, OH and O adsorbates respectively for the model 
that includes the effects of the adsorbate induced rehybridization of the graphene; the parameterization used is given in Table 
HI The solid black lines in (b), (d), (f) and (h) show the calculated conductances of ribbons with H, F, OH and O adsorbates 
respectively for the simpler model that includes only the adsorbate valence orbitals in the EMOs as parameterized in Table lHl 
and thus does not include the effects of the rehybridization of the graphene. Red, green and blue solid lines show the effect on 
the conductance of the individual orbitals of the adsorbed species as indicated. The grey solid line in (a) and (b) shows the 
conductance of a ribbon with interior carbon atom vacancies at p vac = 10 -4 . Ribbon width W = 30 nm; length L = 500 nm. 
Temperature T = 0. t = 2.7 eV. 



of the closest graphene carbon atom, r yH lB ,c 2pz — 1-892 
(see Table Ql]), i s nearly twice as large as the Hamiltonian 
matrix clement t between the 2p z orbitals of adjacent car- 
bon atoms. Adsorption of the H atom also results in a 
modified graphene lattice geometry. In the relaxed geom- 
etry, we find significant lifting of the carbon atom directly 
bound to the H out of the graphene plane by 0.35A. This 
is accompanied by weakening of the C^-p^ — Hamilto- 
nian matrix element between the C atom to which the H 
atom binds and its neighbor C atoms to 0.79t. Since the 
graphene sheet is no longer planar, partial rehybridiza- 
tion from sp 2 to sp occurs near the adsorbed H atom: 
The carbon atom to which the H atom binds can be re- 
garded as forming a bonds with its carbon atom neigh- 
bors and with the H atom. The rehybridization of the 
graphene on adsorption of the H atom and the strong 
coupling between H adatom and the graphene result in 
strong scattering of graphene 7r-band electrons near the 
adsorbed H atom. 



2. Conduction in Ribbons with H Adatoms 

In Fig [5] (a) and (b) we show the calculated conduc- 
tance of a graphene nanoribbon with a concentration 



p H = 10~ 4 of adsorbed H. In Fig [5] (a) the results are 
shown for the model that includes in the adsorbate EMO 
the H Is orbital and the 2s, 2p x and 2p y orbitals of the 
carbon atom to which the H binds and its neighboring 
carbon atoms, the tight-binding parameters used being 
those given in Table |U This model includes the effect 
of the local rehybridization of the graphene from sp 2 to 
sp 3 bonding. For comparison the results for a model in 
which the adsorbate EMO includes only the His orbital 
(with the parameters given in Table lHl of the present pa- 
per) are shown in Fig \5j[b) ■ The conductances of the 
same ribbon with no adsorbate but an equal concentra- 
tion p = 10~ 4 of interior carbon atom vacancies^ and of 
an ideal ribbon with no adsorbate or vacancies are also 
shown. In both models even for this low concentration 
of adsorbed H atoms the conductance of the ribbon is 
strongly suppressed relative to that of the ideal ribbon. 

The calculated conductance with the rehybridization 
of the graphene taken into account (the solid black 
curve in Fig [5] (a)) is strikingly similar both qualitatively 
and quantitatively to that of the ribbon with the same 
concentration of carbon atom vacancies (the solid grey 
curve), whereas the conductance calculated without in- 
cluding rehybridization (the solid black curve in Fig [5] 
(b)) is qualitatively different. 

This is consistent with the idea££ that the rehybridiza- 
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Energy (t) 

FIG. 6: (color online) Comparison of the calculated averaged 
conductances of graphene ribbons with adsorbates (for the 
model that includes adsorbate-induced rehybridization of the 
graphene) and interior carbon atom vacancies at a concen- 
tration p — 10~ 4 in each case. The averaging in each plot 
is over 10 arrangements of the positions of the vacancies or 
adsorbed atoms or molecules. The grey dotted line shows the 
conductance of the ideal ribbon without any defects. 

tion should effectively decouple the carbon atom to which 
the H atom bonds from the graphene n band, in which 
case scattering of graphene 7r band electrons by an ad- 
sorbed H atom would be expected to resemble electron 
scattering by a carbon atom vacancy. However this sim- 
ple picture does not account for the differences between 
graphene with a H adsorbate and graphene with a F or 
OH adsorbate: As will be seen below the conductance 
characteristics of ribbons with F and OH adsorbates with 
rehybridization included differ qualitatively from those of 
ribbons with vacancies and change much less drastically 
than those for H when rehybridization is included in the 
model, although the changes in the graphene geometry 
due to the adsorption of F and OH are very similar to 
and even slightly larger than for H adsorption. 

If the energy scale is broadened (see Fig. |6]) an ap- 
proximately linear increase of the conductance with the 
absolute value of the energy is found with nearly the same 
absolute slope at positive and negative energies for rib- 
bons with adsorbed H. Thus on the larger energy scale 
the conductance of the ribbon is nearly symmetric in the 
Fermi energy in presence of adsorbed H as it is in the 
presence of interior carbon atom vacancies. This behav- 
ior of the conductance is consistent with experiment.— 



3. Role of the Dirac Point Resonance 

A clearer understanding of the role of the graphene 
rehybridization in electron transport can be gained by 



comparing the the conductance plots in Fig [5] (a) (where 
the rehybridization is included in the model) and Fig 
[5] (b) (where it is not) with the calculated properties of 
the Dirac point scattering resonances associated with ad- 
sorbed H in the same two models as were used for these 
transport calculations. The corresponding Dirac point 
resonances are shown in Fig. [2] and [3] respectively. Ignor- 
ing the mesoscopic conductance fluctuations in Fig. [SJ 
it is apparent that including the rehybridization in the 
transport model results in a shift in the energy of the 
main conductance minimum from ~ — O.li in FigJSJb) 
to ~ O.Ot (the Dirac point of graphene) in FigJSJa). This 
matches reasonably well the shift in the energy of the hy- 
drogen Dirac point scattering resonance from — 0.138i to 
— 0.0026i that was found upon inclusion of rehybridiza- 
tion in the theory of the H Dirac point resonance in Sec- 
tion [VT] Thus the locations in energy of the conductance 
minima for the two models agree quite well with the en- 
ergies at which the Dirac point resonances of H occur 
in the models. This is reasonable since strong resonant 
electron scattering near a particular energy is expected 
to suppress electron transport near that energy. As will 
be seen below there is similar agreement between the 
Dirac point resonance energies and the energies at which 
the nanoribbon conductance minima occur for the other 
adsorbates that we study in this article. 

A notable feature of the conductance calculated for the 
nanoribbons with adsorbed hydrogen in the model that 
does not include the graphene rehybridization (FigJSJb)) 
is that in this case the conductance is affected only 
weakly by the presence of adsorbate when only the lowest 
nanoribbon subband is populated with electrons. That 
is, the conductance of the ribbon with the H adsorbate 
in FigJSJb) near the Dirac point is very close to 2e 2 /h, 
the conductance of the pristine ribbon in the same en- 
ergy range. This remarkable robustness of transport in 
the first subband is due to that subband's unique scat- 
tering properties arising from the nature of the uncon- 
fined electron wave function in graphenei 25 ' 35 Upon in- 
clusion of the graphene rehybridization in our model we 
found in Section IVII that in addition to the Dirac point 
scattering resonance for the H adsorbate shifting to an 
energy very close to the graphene Dirac point the reso- 
nance also becomes stronger by approximately two orders 
of magnitude. This is sufficient to override the relative 
robustness of transport in the first graphene nanoribbon 
subband against scattering so that, unlike in FigGJb), 
the conductance close to the Dirac point (zero energy) 
in FigJSJa) for the ribbon with the H adsorbate is very 
strongly suppressed relative to that for the ideal ribbon. 



4- Adsorbed H Atoms and the Density of States 

The presence of the adsorbate leads to electron local- 
ization at the adsorbed atom or molecule as well as at the 
nearby carbon atoms. Fig. [7] (a) shows the local density 
of states (LDOS) for a ribbon with H adatoms averaged 
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over atoms at different locations. Electron localization is 
most pronounced for energies near the Dirac point where 
the conductance is suppressed most strongly due to res- 
onant backscattering of electrons. Localization is much 
stronger at the H adatom (including its EMOs) than at 
the 2p z orbital of the carbon atom to which the H binds. 
The LDOS for that carbon 2p z orbital is smaller by three 
orders of magnitude than for the H atom. For this rea- 
son it is indistinguishable from the horizontal axis in Fig. 
[7J (a). The carbon atoms adjacent to that carbon atom 
belong to the other graphene sublattice and exhibit an 
LDOS as large as that at the H. Strong electron local- 
ization is known to occur in graphene near carbon atom 
vacancies . 33 ' 38 ' 59 The lower left and middle panels in Fig. 
[Jj show the LDOS for two representative graphene rib- 
bons with interior vacancies and H adatoms, with these 
defects located at corresponding sites in the two ribbons. 
For the chosen energy — 0.03i electrons localize strongly 
near defects of both types. The C3 point symmetry of the 
graphene lattice about the defect sites is clearly visible in 
LDOS in the vicinity of each defect. We find the LDOS 
to decay according to the power law 1/r (not shown in 
Fig. [7]) for both defect types although individual defects 
may exhibit differing LDOS amplitudes. This power law 
decay is consistent with the results of other studiesi 67 ' 103 



B. Ribbons with F Adatoms 

The fluorine adatom bonds covalently to a carbon atom 
of the graphene lattice and has two valence orbitals (2s 
and 2p z ) that scatter graphene 7r electrons. In the sim- 
plest tight-binding model that ignores rehybridization of 
the graphene, these two orbitals represent two indepen- 
dent channels for electron scattering and their contribu- 
tions to the scattering are additive (for a single adatom). 
However, the 2s and 2p z orbitals of the F do not con- 
tribute equally: The 2s orbital couples more strongly to 
the carbon 2p z orbital since it has the larger value of 
\la,c 2p , I as can be seen in Tabic HT1 It therefore has the 
most influence on the electron transport in the graphene 
ribbon. The coupling of the 2p z orbital is much weaker 
and has a much smaller effect on the conductance; see 
Fig. ETd). 

The Dirac point resonance energy for F on graphene 
was found in Section fVTl to be — 0.1 36i in the model that 
includes rehybridization and — 0.222t in the model that 
does not. Thus while the effect of rehybridization of the 
graphene on the Dirac point resonance of F is important 
it is not as drastic as in the case of H. However, as in the 
case of adsorbed H, these Dirac point resonance energies 
for F agree well with the energies at which the calcu- 
lated conductances of ribbons with adsorbed F are most 
strongly suppressed in both models as can be seen in Fig. 
HHc) and (d) and Fig. © 

In contrast to the case of the H adsorbate, the con- 
ductance is affected only weakly by the presence of the 
F adsorbate when only the lowest nanoribbon subband 
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FIG. 7: (color online) Local density of states (LDOS) aver- 
aged over 50 adatoms located at different places in the ribbon 
vs. electron energy for a ribbon with adsorbed (a) H and (b) 
F atoms. Adsorbate-induced rehybridization of the graphene 
is included. The LDOS is shown in red for the 2p z orbital of 
the C atom to which the adsorbed atom bonds and in green 
for the 2p z orbital of a C atom neighbor of that C atom. 
The total LDOS on the EMOs associated with the adsorbed 
atom (including the contributions of the non-7r orbitals of the 
nearby C atoms) is shown in blue. The dotted and solid grey 
lines show the conductances of the ideal ribbon and the rib- 
bon with H or F adatoms, respectively. Arrows mark the 
energy chosen for computation of the LDOS in the bottom 
panels. The latter are for 15 x 30 nm 2 fragments of ribbons 
with interior vacancy defects (left), H adatoms (middle) and 
F adatoms (right) located at the same places. In the mid- 
dle and right panels the LDOS for the adsorbate EMOs is 
combined with the graphene LDOS. For the chosen energy 
(— 0.03t) the missing carbon atoms give rise to very similar 
electron localization close to the defects as the presence of the 
H adatoms while the LDOS for F is noticeably different. 



is populated with electrons for both models at the adsor- 
bate concentrations in Fig. [Sfc) and (d). This is because 
(unlike for H) the Dirac point resonances for F are off- 
set significantly in energy from the Dirac point in both 
models. 
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The calculated LDOS for the F adsorbate for the 
model that includes the rehybridization of the graphene 
is shown in FigJTJb) for the EOMs asociated with the F 
adsorbate in blue and for the 2p z orbitals of the C atom 
to which the F binds (red) and its C neighbor (green). 
In each case the LDOS is strongest around the energy of 
the F Dirac point resonance. As for H the LDOS for F 
is weaker on the 2p z orbitals of the C atom to which the 
F binds than either on the neighboring C atoms or on 
the EMOs of the F adsorbate. The spatial map of the 
LDOS for F shown in the lower right panel of FigJT] is 
qualitatively similar to that for H and for vacancies, but 
the particular defect sites showing the strongest LDOS 
at a given energy are in some cases different for the F. 

Notice also that in Fig[7fb) there are peaks in the 
DOS at positive energies near subband edge energies that 
match dips in the conductance that is also plotted for 
comparison. These conductance dips are due to enhanced 
electron backscattering due to the greater availability of 
final states for the scattering process at energies with the 
higher density of subband states as is discussed in Ref. 



C. Ribbons with Adsorbed OH Groups 

The OH molecule is another monovalent adsorbate 
with many properties similar to F adatoms. In the re- 
laxed geometry the O bonds covalently to a C atom (top 
site bonding) and the OH chain stands upright relative to 
the graphene plane, features found also in Refs. Islsi 
As is seen in Fig. [5jf ) , the 2s and 2p z O orbitals af- 
fect electron conduction through the graphene ribbons 
similarly to 2s and 2p z orbitals of F (see Fig. [5jd)) 
although for O the orbitals scatter electrons somewhat 
more strongly. The H Is orbital in OH molecule has 
much less influence because of its large distance from the 
C and small overlap with the 2p z C orbital. The to- 
tal effect of OH adsorbed molecules on the conductance 
through the graphene ribbon is qualitatively and quanti- 
tatively similar to that of F; see also Fig. [6] 

The Dirac point resonance energy for OH on graphene 
was found in Section Ell to be — 0.089i in the model that 
includes rehybridization and — 0.194i in the model that 
does not. These numbers again agree reasonably well 
with the energies at which the conductances of ribbons 
with adsorbed OH is strongly suppressed. The Dirac 
point resonances for OH occur at energies intermediate 
between those for H and F in each model and the same 
is true of the energies at which the strongest suppression 
of the calculated conductance is seen in Fig. [5] and [6] 

As in the case of the F adsorbate, the conductance is 
affected only weakly by the presence of the OH adsorbate 
when only the lowest nanoribbon subband is populated 
with electrons for both models at the adsorbate concen- 
trations in Fig. [Sfe) and (f), for the same reasons. 



D. Ribbons with O Adatoms 

Oxygen is a bivalent adsorbate that binds simultane- 
ously to two neighboring carbon atoms (a bridge site). 
This leads to substantial rehybridization of the bonding 
associated with these two carbons that belong to differ- 
ent graphene sublattices. This, in turn, leads to strong 
electron scattering and suppression of electron conduc- 
tion through the ribbon. 

The calculated conductance vs. energy characteris- 
tics for ribbons with adsorbed O are shown in Fig. [5jg) 
and [5] for the model that includes rehybridization of the 
graphene and in Fig. [5{h) for the model that does not. 
The O orbitals affect conduction differently: The O 2s 
orbital produces strong suppression of the conductance 
at positive energies, but the 2p z orbital suppresses the 
conductance much less and this occurs mainly at nega- 
tive energies. The effect of the O 2p x orbital is weaker 
still. The resulting low conductance at positive energies 
is a consequence in part of the O adsorbate binding to 
two carbon atoms belonging to two graphene sublattices 
unlike the F and OH that bind to a single C atom and 
exhibit low conductance at negative energies. We note 
however that in general the sign of the energy at which 
low conductance due to adsorbate scattering occurs de- 
pends not only on the number of C atoms to which the 
adsorbate binds but also on the values of the model tight 
binding parameters e a and j a j. For example, if the ab- 
solute values of 7 aj - were substantially smaller than our 
estimate the low conductance region for O would be at 
negative energies relative to the Dirac point of the ribbon. 

As for the other adsorbates discussed above both the 
sign and magnitude of the energy at which the conduc- 
tance of the ribbon with the O adsorbate is suppressed 
agree well with the energies at which the Dirac point 
resonances for the O adsorbate were found to occur in 
Section EH i.e., 0.112* and 0.090i in the model that in- 
cludes the graphene rehybridization and 0.046£ for that 
which does not. 

Another striking feature of Fig. [5] is the prominent 
conductance maximum for the ribbon with the O adsor- 
bate at energies near 0.55i. The energy at which the 
conductance maximum occurs coincides with the energy 
of the antiresonance (i.e., deep minimum) of the norm 
of the T-matrix (that describes scattering of graphene 
electrons due to an adsorbed O atom) that was found in 
Section IVII in the model that includes the O adsorbate- 
induced rehybridization of the graphene. 

Since weak electron scattering is normally associated 
with high conductance this match between the conduc- 
tance maximum and the T-matrix antiresonance is in- 
tuitively reasonable. It demonstrates once again the 
close relationship between the T-matrix theory of Sec- 
tions [IV] and EH and transport in graphene ribbons with 
chemisorbed species that has already been illustrated by 
the agreement found above between the electron energies 
at which T-matrix resonances and strong suppression of 
the ribbon conductances occur. 
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E. Conductance Asymmetry Relative to the Dirac 

Point of Graphene 

In Section WT\ we found the Dirac point resonances for 

F, OH and to be offset in energy from the Dirac point of 
graphene and in Sections IIX B 1 IPX~Cl and ITX D I we showed 
the electron Fermi energies at which the lowest conduc- 
tances of ribbons with these adsorbates occur to be offset 
from the Dirac point of graphene similarly. 

Such asymmetric conduction relative to the Dirac 
point has been discussed in Ref. [H] for the case of in- 
finite two-dimensional graphene. Suppression of the con- 
ductivity on one side of the Dirac point with only a weak 
effect on conduction on the other was explained by a 
local energy- dependent scattering potential due to the 
adsorbate^ Over a range of energies the conductivity 
was found to be small while rising linearly outside that 
range. Qualitatively similar results for 2D graphene were 
subsequently reported by othersj 70 ' 71 A transport gap 
around the impurity resonance energy had also been pre- 
dicted earlier for a generic model of infinite 2D graphene 
with point defects 60 although this prediction was not sup- 
ported by transport calculations^ 

This behavior resembles our results for the conduc- 
tance in the graphene ribbons with F, OH and O ad- 
sorbates but with a qualitative difference: As will be 
discussed in Section IIX Ft at moderately low temper- 
atures the graphene ribbons exhibit quantized conduc- 
tance steps superposed on the otherwise linearly rising 
conductance due to enhanced electron backscattering at 
the edges of the subbands of the ribbon. 



F. Conductance Quantization 

In the experimental study of Lin et al£ conductance 
quantization in the form of conductance steps of equal 
height was observed in graphene nanoribbon samples 
with conductances much smaller than 2e 2 /h in a range of 
moderately low temperatures as a gate voltage applied to 
the sample was varied. In a previous paper— we showed 
theoretically that conductance quantization of this kind 
should occur in graphene nanoribbons with interior car- 
bon atom vacancies even if comparable amounts of other 
defects such as edge disorder and long range potentials 
due to charged defects are also present. However, be- 
cause of the sample preparation techniques used by Lin 
et alj^ adsorbed H may well have been present in their 
samples. 

As we have already noted in Section IIX A 21 (and is 
demonstrated very clearly in Fig [5] (a)) the calculated 
zero temperature conductance of graphene ribbons with 
adsorbed H atoms is very similar both qualitatively and 
quantitatively to that of the ribbon with the same con- 
centration of carbon atom vacancies. In particular, the 
properties of the conductance characteristics of ribbons 
with vacancies that have been shown4°- to give rise to 
conductance quantization of the kind observed by Lin et 
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FIG. 8: (Color online) The calculated conductances of ideal 
graphene ribbons with no adsorbate (solid grey) with H 
adatoms (solid black) and with F adatoms (dashed red). The 
temperature is 0, 80 and 300 K in (a), (b) and (c) respectively. 
Adatom concentrations are p = 10 -4 as in Fig[5] Adsorbate 
induced graphene rehybridization is included in the model. 
Ribbon width W = 30 nm; length L = 500 nm. t = 2.7 eV. 



alj£ are also exhibited by the conductance characteris- 
tic of the ribbon with H adatoms in Fig [5] (a) . These 
properties are 

(i) Pronounced sample-specific conductance fluctua- 
tions that are manifestation of quantum interference 4°- 

(ii) If the conductance fluctuations are ignored, the 
conductance is seen to scale down uniformly overall due 
to scattering by the adsorbate, i.e., in a similar way for 
all subbands. 

(iii) The conductance shows a pronounced dip when- 
ever a new subband becomes available for electron prop- 
agation. As is discussed in Ref. 40, this is because of 
enhanced electron backscattering by the defects at sub- 
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band edges. 

Because of these three properties we expect graphene 
nanoribbons with adsorbed H to exhibit equally spaced 
conductance steps similar to those observed experimen- 
tally by Lin et air- for the same reasons and under the 
similar conditions (discussed in Ref. H3) as do graphene 
nanoribbons with carbon atom vacancies. That is, con- 
ductance steps of equal height should be observed even 
in samples with sufficiently high adsorbate concentra- 
tions for the ribbons to have conductances much smaller 
than 2e 2 /h, the conductance steps should break up into 
random conductance fluctuations as the temperature ap- 
proaches zero Kelvin, and the conductance steps should 
become completely smeared out by thermal broadening 
at temperatures substantially larger than the subband 
spacing of the ribbons. This is indeed seen in Fig. [8] 
where we show the calculated conductance of a ribbon 
with H adatoms at different temperatures: Regular con- 
ductance features are seen at 80K in Fig. [5] (b). They 
break up into universal conductance fluctuations at OK 
in Fig. [5] (a) and are completely smeared out thermally 
at 300K in Fig. H](c). 

Thus electron scattering due to a low concentration 
of H atoms adsorbed on the graphene nanoribbons in 
our model that takes into account the adsorbate-induced 
rehybridization of the graphene provides an alternative 
and equally satisfactory explanation of the conductance 
quantization observed by Lin et ali^ to thats* provided 
by electron scattering by carbon atom vacancies in the 
interior of the ribbon. 

It is also possible that adsorbed F, OH and O were 
present in the samples of Lin et al£ The three properties 
(i), (ii) and (hi) are also shared by ribbons with F, OH 
and O adsorbates for positive, positive and negative en- 
ergies respectively, as can be seen in Fig. [S](c), (e) and 
(g). (Notice also the density of states maxima at sub- 
band edges that are responsible 4 ^ for enhanced electron 
backscattering and hence for a conductance dip when- 
ever the electron Fermi level crosses a subband edge (i.e., 
property (iii) above) that are clearly visible at positive 
energies in Fig. EJb)). Therefore, conductance quanti- 
zation of the kind observed by Lin et al*2, should also 
occur under appropriate conditions for some ranges of 
the gate voltage in graphene ribbons with adsorbed F, 
OH and O. Theoretical results for H and F adatoms at 
various temperatures are shown in Fig. [5J Note that 
the nearly perfect transmission of electrons in the first 
subband through the ribbons with F, OH and O seen in 
Fig. [5] and [7] is specific to armchair ribbons of partic- 
ular widths^! that are metallic and does not occur for 
insulating ribbons. However, unlike for adsorbed H, for 
adsorbed F, OH and O, the Dirac point has a distinc- 
tive signature in the calculated conductance that does 
not coincide with the conductance minimum that is due 
to the scattering resonance(s) associated with the adsor- 
bate. Such a signature is not evident in the experimental 
data of Lin et alA 

In the experimental data of Lin et al<£ there is an off- 



set of about 3 Volts between the gate voltage at which 
the conductance minimum occurs and zero gate voltage. 
The offsets between the energies at which the conduc- 
tance minima occur in our calculations and the Dirac 
point energy for F, OH and O are all smaller than 0.3 
eV, the largest being for the F adsorbate. However, we 
expect other mechanisms such as charged impurities in 
the insulating spacer between the ribbon and gate elec- 
trode, contact potentials between the spacer and gate 
and/or between the spacer and ribbon to contribute sig- 
nificantly to the the experimentally observed gate voltage 
at which the conductance minimum occurs. 



G. Adsorbate Induced Renormalization of the 
Dirac Point and Subband Edge Energies 

We note that while for F, OH and O the conductance 
minimum is displaced in energy from the Dirac point of 
the pristine ribbon, this should not be interpreted as a 
shift of the Dirac point of the ribbon due to interaction 
with the adsorbatejiS 4 - it arises almost entirely from sup- 
pression of the conductance by enhanced electron scat- 
tering near the conductance minimum. 

The actual shift of the Dirac point energy (and sub- 
band edge energies) of the ribbon due to interaction of 
the ribbon with the adsorbate is much smaller and can 
be estimated perturbatively as follows: If we consider 
the coupling term H c — J2 a j Jaj (d^aj + /i.e.) between 
the adsorbate and the 7r states the ribbon in the tight 
binding Hamiltonian Eq. (fT5j) as a perturbation then in 
second order perturbation theory this coupling implies a 
shift Afc s in the energy of the subband state \ks) given 
by 

A fcs - J2\^\ H c\ks)\ 2 /(e ks - e a ) (19) 

where efe s is the unperturbed energy of state \ks). For 
states \ks) that are much closer in energy to the Dirac 
point of the ribbon than are the orbital energies e Q of 
the adsorbate we can approximate tk s — £q ~ — in the 
denominator of Eq. fTS] For OH and F adsorbates that 
bind to the graphene over a single C atom and assuming 
that they are randomly distributed over the ribbon we 
can on average approximate \(ip a j\H c \ks)\ 2 by |7aj| 2 /^ 
where N is the number of carbon atoms in the ribbon. 
With these approximations Eq. fT9l becomes 

A fes ~ -p^lTajfAa (20) 
a 

where p is the concentration of the adsorbate and the 
sum is over the extended molecular orbitals of a single 
adsorbed atom or molecule. Inserting p — 10~ 4 and the 
values of j a j and e a from Table Q] in Eq. [201 we find 
A fcs ~ 4.4 x IQ~H and 6.3 x 10~ 4 t for F and OH re- 
spectively These shifts are very small justifying our use 
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of perturbation theory and demonstrating that the shift 
of the Dirac point energy due to this mechanism is very 
small for the adsorbate concentrations considered in this 
article. 

Despite its small size the energy shift Afc s has a clear 
signature in the results of our transport calculations: As 
we have discussed in Ref. and in the preceding subsec- 
tions of Section llXl of the present article, the conductance 
of a ribbon with disorder (including that due to ran- 
domly located adsorbed moieties) shows a pronounced 
dip whenever a new subband becomes available for elec- 
tron propagation due to enhanced electron backscatter- 
ing by the defects at subband edge energies. Given that 
the interaction between the ribbon and adsorbate renor- 
malizes the subband edge energies by the amount Afc s 
it is to be expected that the energies at which the con- 
ductance dips due to scattering by an adsorbate occur 
should be offset from the subband edge energies of the 
pristine ribbon also by approximately Afc s . This is in- 
deed what we find. For example, for F adatoms on a 
ribbon at a concentration p = 4 x 10 -4 Eq. [501 yields 
Afc s =~ 1.8 x 10~ 3 £ which agrees reasonably well with 
the displacement D ~ 1.5 x 10~ 3 t of the center of the 
conductance dip from the second subband edge of the 
pristine graphene ribbon that is seen for the p = 4x 10~ 4 
F adsorbate in Fig. EJc). 



H. The Transport Gaps and Their Dependence on 
Adatom Concentration 

Figure [9] shows the conductances of graphene ribbons 
with different concentration of H and F adatoms. As 
the adsorbate concentration increases the conductance 
decreases and a wider transport gap in which the ribbon 
is effectively an insulator opens centered near the Dirac 
point for H and at negative energies for adsorbed F. The 
width A g of the transport gap (that we define arbitrarily 
as the energy range where the ribbon conductance is less 
than 0.9 x 2e 2 /h) grows linearly with the adatom con- 
centration at low concentrations. It worth noting that 
A ff does not depend on the ribbon width, which rules 
out adsorbate impurities as the principal source of the 
transport gap in the experiment in Ref. QJ. 

Transport gaps gaps have been discussed previously in 
Ref. 1 1051 for graphene with uncompensated vacancies. A 
A g y/p dependence was found over a wide range of 
defect concentrations < p < 0.2. Note, however, that 
for vacancies randomly distributed over the two sublat- 
tices no transport gap was found in Ref. Il05l Transport 
gaps have also been predicted for graphene ribbons sub- 
stitutionally doped with boron,— however, whether the 
predicted transport gaps are related to Dirac point reso- 
nances was not discussed 4^ 




P Energy (t) 




Energy (t) 

FIG. 9: The averaged conductances of graphene ribbons with 
different concentrations p of adsorbed atoms of (b) H and (c) 
F averaged over different locations of the adsorbed atoms. 
The results for p = 0.5 x 1(T 4 , 1 x 1(T 4 , 2 x 1(T 4 and 4 x 1(T 4 
are shown in green, blue, red and black, respectively. The av- 
eraging in each plot is over 10 arrangements of the positions 
of the adsorbed atoms. The transport gap A 9 is estimated as 
the energy interval where G < 0.9 x 2e 2 /h. (a) A 9 increases 
approximately linearly with p at low p. D ~ 1.5 x 10 -3 t in (c) 
is the offset between the second subband edge for the ideal rib- 
bon with no adsorbate and the center of the conductance dip 
due to enhanced electron backscattering when the (renormal- 
ized) edge of the second ribbon subband crosses the electron 
Fermi level in the ribbon with a concentration p = 4 x 10 -4 
of adsorbed F atoms; see the last paragraph of Section IIX Gl 
Note that p — 10 -4 corresponds to 3.8 x 10 15 adsorbed atoms 
per square meter, t — 2.7 eV. 



X. DISCUSSION 

In this article we have formulated a tight binding the- 
ory of the Dirac point resonances due to adsorbed atoms 
and molecules on graphene based on the standard tight 
binding model of the graphene 7r-band electronic struc- 
ture and the extended Huckel model of the adsorbate 
and the adsorbate-induced local sp 3 rehybridization of 
the graphene. We generalized previous theories of the ef- 
fective Hamiltonians of graphene with impurities to the 
case of adsorbate species with multiple extended molec- 
ular orbitals and bonding to more than one graphene 
carbon atom, and obtained accurate analytic expressions 
for the Green's function matrix elements that enter the 
T-matrix theory of Dirac point resonances. This gener- 
alization makes the extended Huckel model (and poten- 
tially other tight binding models as well) into a power- 
ful tool for studying the Dirac point resonances induced 
by many different adsorbatcs on graphene. Furthermore 
this theory makes it practical to carry out sophisticated 
electronic quantum transport calculations for graphene 
nanoribbons tens of nanometers wide (such as are being 
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realized in present day experimental studies) with ad- 
sorbates covalently bonded to the ribbon. We applied 
the above theoretical approach to H, F, OH and O ad- 
sorbates on graphene whose relaxed geometries we calcu- 
lated with ab initio density functional theory. For each of 
these adsorbates we found a strong scattering resonance 
near the Dirac point of graphene, the strongest by far 
being for the H adsorbate. Treating the valence orbitals 
of the adsorbed species and the 2s, 2p x and 2p y valence 
orbitals of the nearby carbon atoms theoretically in a uni- 
fied way was necessary in order to obtain reliable results. 
We also extracted from these calculations a minimal set 
of tight binding parameters that make it possible to effi- 
ciently model adsorbate-induced electron scattering and 
its effect on electron transport in graphene and graphene 
nanostructures. 

In particular, the minimal tight binding models that 
we developed make it possible to model the effect of ad- 
sorbates on transport in graphene nanoribbons tens of 
nanometers wide and hundreds of nanometers long that 
are at present being realized experimentally. We have 
presented realistic electronic quantum transport calcula- 
tions for such nanoribbons with adsorbed H, F, OH and 
O. As well as the carbon 7r band electronic structure of 
the graphene nanoribbons our theory includes the effects 
of the local partial rehybridization of the graphene ribbon 
from the sp 2 to sp electronic structure that occurs when 
H, F, OH or O bonds covalently to the ribbon. This is 
necessary in order for the model to describe correctly the 
scattering resonances that are induced in the graphene 
ribbons near the Dirac point by the presence of these 
adsorbates. We find that these Dirac point resonances 
play a dominant role in quantum transport in ribbons 
with these adsorbates: Even at low adsorbate concentra- 
tions of 10~ 4 adsorbed atoms or molecules per carbon 
atom, in the ribbons that we study the conductance of 
the ribbon is strongly suppressed and a transport gap is 
formed for electron Fermi energies in the vicinity of the 
energy of the resonance. For the H adsorbate this trans- 
port gap is centered very close to the Dirac point energy 
of the ribbon (as it is for ribbons with interior carbon 



atom vacancies) while for F and OH it is centered below 
the Dirac point and for O it is centered above the Dirac 
point. These predictions can be tested experimentally 
by tuning the position of the Fermi level in the ribbon 
relative to its Dirac point by means of a variable applied 
gate voltage. For each of these adsorbed species we find 
a pronounced dip in the low temperature conductance 
of the ribbon when the electron Fermi level crosses the 
edge of an electronic subband of the ribbon due to en- 
hanced electron backscattering, and the conductance to 
be suppressed equally on average in every subband. This 
implies that graphene nanoribbons with H, F, OH and O 
adsorbates and conductances even a few orders of magni- 
tude smaller than 2e 2 /h should (for appropriate ranges of 
a back gate voltage) exhibit equally spaced conductance 
steps at moderately low tempertaures similar to those 
that have been observed by Lin et al.— experimentally 
and that have recently been predicted theoretically 4 ^ for 
graphene nanoribbons with interior carbon atom vacan- 
cies. 

Experiments testing these predictions by observing the 
Dirac point resonances in lateral transport through well 
characterized graphene nanoribbons with H, F, OH and 
O adsorbates intentionally deposited at known concen- 
trations would be of interest. Vertical transport measure- 
ments directed at detecting the Dirac point resonances 
and measuring the energies at which they occur more di- 
rectly with the help of scanning tunneling spectroscopy of 
atoms and molecules adsorbed on graphene and graphene 
ribbons would also be interesting, especially in view of 
the often conflicting predictions of the values of some of 
these energies that have been made by various theory 
groups. 
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